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We study the effect of electron-electron interaction and spin on electronic and transport properties 
of gated graphene nanoribbons (GNRs) in a perpendicular magnetic field in the regime of the lowest 
Landau level (LL). The electron-electron interaction is taken into account using the Hartree and 
Hubbard approximations, and the conductance of GNRs is calculated on the basis of the recursive 
Greens function technique within the Landauer formalism. We demonstrate that, in comparison to 
the one-electron picture, electron-electron interaction leads to the drastic changes in the dispersion 
relation and structure of propagating states in the regime of the lowest LL showing a formation of 
the compressible strip and opening of additional conductive channels in the middle of the ribbon. 
We show that the latter are very sensitive to disorder and get scattered even if the concentration of 
disorder is moderate. In contrast, the edge states transport is very robust and can not be suppressed 
even in the presence of a strong spin-flipping. 

PACS numbers: 72.80.Vp, 73.22.Pr, 73.43.Cd, 73.43. Qt 



I. INTRODUCTION 

Due to Dirac spectrum graphene possesses a num- 
ber of unique electronic and transport properties^. In 
a magnetic field B perpendicular to a graphene layer 
the spectrum is modified into a series of Landau levels 
(LL) with the energies E = ±hiOc\/N, where the cy- 
clotron frequency ujc = vpy^2eB/h, vp ~ 10* cm/s is 
the Fermi velocity and = 0, ±1, ±2, . . . ^ One of inter- 
esting features of this spectrum is a presence of O'th LL 
with E = 0, which is equally shared between electrons 
and holes. In contrast to conventional two-dimensional 
electron gas systems the position of this level does not 
depend on a value of the magnetic field. An experimen- 
tal manifestation of these unusual series is the anomalous 
quantum Hall effect with the Hall conductivity given by 
o'xy = + 1/2), where the factor 4 comes from the 

spin and valley degeneracy^il. In a strong magnetic field 
the four- fold degeneracy of the lowest (i.e. O'th) LL is 
lifted leading to the insulating behavior characterized by 
the presence of additional plateaux in the Hall conductiv- 
ity f^xy accompanied by the increase of the longitudinal 
resistivity Pxxi^ — The large energy gaps identified experi- 
mentally suggests that electro-electron interactions along 
with the Zeeman spin-splitting play an important role in 
understanding of this phenomena. 

Most of the experiments reported so far are performed 
in the four-terminal geometry which allows a direct mea- 
surement of the components of the resistivity and con- 
ductivity tensors^"—. However in some cases the only 
possible experimental setup is the two-terminal measure- 
ments, even though the results obtained in this geome- 
try are not as straightforward to interpret as the four- 
terminal ones^. A difference between the longitudinal 
and the Hall resistances is not clearly defined in such the 
measurements. The unexpected consequence of analysis 
of the two-terminal geometry is that a measured two- 
terminal resistance between the right and left contacts, 
i?2t, corresponds, in the four-terminal setup, to the Hall 



resistance Rh rather than to the longitudinal one Rl- 
This stems from the fact that the voltage drop between 
the sample edges, V, equals to the difference of the chem- 
ical potentials in the contacts eV = /-t/j — hl J^^ii 

An explanation of the quantum Hall effect is often 
based on the picture of edge states. In a strong enough 
magnetic field the right and left propagating states are lo- 
calized at different edges of the system due to the Lorentz 
force. It leads to exponentially small overlap between 
counterpropagating states which, in turn, greatly sup- 
presses a possibility of backscattering and is eventually 
manifested by developing pronounced plateaux in the 
conductance. In graphene another transport regime can 
be realizedi^. Coupling of spin to magnetic field due to 
the Zeeman effect leads to a splitting of the lowest LLs 
into two sublevels. For the neutral graphene the chem- 
ical potential coincides with the zero energy and lies in 
the center of the spin-gap, see Fig. [3] (d) below. At a 
given Fermi energy there are two forward propagating 
states of opposite spins located near different edges. The 
most prominent feature of this transport regime is that 
forward and backward propagating states of the different 
spins are not spatially separated. If a scattering event 
leads to a spin fiipping, the overlap between the counter- 
propagating states of different spins might be sufficient 
to cause the backscattering. This backscattering of the 
spin-polarized states is used for an explanation of the 
non-zero behavior of pxx near the Dirac pointi^. 

It should be stressed that the structure of edge states 
derived within a one-electron approximation and shown 
on Fig.([2|) (d) corresponds to the case of neutral graphene 
when the electron density is zero. It is often used for 
understanding and explanation of a variety of experi- 
ments on magnetotransportji^. The conditions for ex- 
istence of the edge states in GNRs was discussed in de- 
tails in Ref.fi^]. Also, self-consistent Hartee-Fock calcu- 
lations were performed for GNR in magnetic fieldi^. It 
was shown than the favored ground states in the neu- 
tral graphene are charge- or spin-density waves. Neither 
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of these orderings support the edge states. However, in 
experiments one usuahy needs to tune a charge concen- 
tration by means of the gate voltage. It was shown be- 
fore that electron-electron interaction between induced 
charges can dramatically change the energy dispersion 
and the structure of propagating st&tesi^r—. The mag- 
netosubband structure and the character of propagating 
states in the lowest LL for the case of nonzero induced 
charge density remain largely unexplored and represent 
the main focus of the present study. We investigate 
the effect of the electron-electron interaction, spin and 
global gate electrostatics on magnetotransport proper- 
ties of GNR in the vicinity of the charge neutrality point. 
Based on the calculated self-consistent magnetosubband 
structure we compute and analyze the two-terminal con- 
ductance of realistic GNRs focussing on the effect of dis- 
order, surface warping and spin flipping on the robustness 
of the edge state transmission. 



II. MODEL 
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FIG. 1: (Color online). A schematic layout illustrating a sys- 
tem consisting of a graphene nanoribbon of the width w = 20 
nm which is located on an insulating substrate at the distance 
d = 300 nm apart from a metallic back gate. The whole sys- 
tem is subjected to a perpendicular magnetic field B — 150 
T. 



We consider a system consisting of a GNR of the width 
w = 20 nm and the length I = 85 nm, located on an 
insulating substrate (e^ = 3.9) and connected to ideal 
graphene leads of the same width, see Fig. \T\ A metallic 
back gate is used to tunc the Fermi energy and the charge 
concentration in the GNR. The whole system is subjected 
to a uniform magnetic field, B = 150 T, perpendicular 
to the graphene plane. Such the high value of the mag- 
netic field allows us to use relatively narrow GNR (in 
order to reduce computational time) and choose the ra- 
tio w/Zb sa 9.5 in an accordance with typical experiments 
(with Ib = y^h/eB being the magnetic length). In order 
to model magnetotrasport in a GNR we use the tight- 
binding Hubbard Hamiltonian in the mean-field approx- 
imation H = -\- which is shown to describe carbon 
electron systems in good agreement with first-principles 
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where a, a' describe two opposite spin states f: 4-5 
the summation runs over all sites of the graphene lat- 
tice, A includes the nearest neighbors only, tr,r-(-A = 
to exp(i27r(/)r,r-i-A/'/'o) with = 2.11 eV, c/iq = hje be- 
ing the magnetic flux quantum which is calculated using 
the Pierel's substitution 



Pr,r+A 



(2) 



with A being the vector potential. (In our calculations 
we use the Landau gauge, A = (— i3y,0)). The first 
term in the Hamiltonian ([T]) describes hopping of elec- 
trons through the graphene lattice. The second term 
describes the spin-splitting in the magnetic field B due 
to the Zeeman effect. Both these terms correspond to the 
non-interacting part of the Hamiltonian. The interaction 
between electrons of the opposite spins is described by 
the Hubbard term with U = 1.33t(r^i2i, and V{r) is a 
Hartree term describing the long-range Coulomb inter- 
action between induced charges in the GNRiiSi^S, 
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- r'|2 + 4^2 / 
(3) 

where d is the distance between the GNR and the gate, 
and the second term describes a contribution from the 
mirror charges. In order to calculate the electron number 
at site r, Eq.([T]) is solved using the recursive Green's 
function technique^ 
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{rirrr = 



[G„{r,r,E)]fFDiE,EF)dE, (4) 



where G'o-(r, r, E) is the Green function in the real space 
representation of an electron of the spin a residing on 
the site r and Ef = eVg is the Fermi energy which po- 
sition is adjusted by the gate voltage. Equations 
are solved self-consistently in order to compute the band 
structure, the charge density and the potential distribu- 
tion. For a given potential distribution we compute the 
conductance using the Landauer formula 



GiEF,B)^Go J T{E,B) 



dfFpjE-EF) 
dE 



dE, (5) 



where T{E, B) is a total transmission coefficient, Gq = 
/h is the conductance quantum, and Jfd is the Fcrmi- 
Dirac distribution. 
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FIG. 2: (Color online), (a)-(c) Charge density distributions, and (d)-(f) band structure of a graphene nanoribbon in a 
perpendicular magnetic field B — fSOT. (a),(d) one-electron approximation, eVg/to = 0; (b),(e) self-consistent calculations for 
eVg/to — 0.1 and (c),(f) eVg/to = 0.5. Red and blue curves marked by t and 4- correspond to the charge densities of spin-up 
and spin-down electrons respectively, nf{y),ni{y). Charge densities are averaged over three successive sites. Green curves 
correspond to the total density distribution n{y) = n^{y) +'ni{y). Dashed lines define the Fermi energy position. Yellow fields 
correspond to the energy interval [—2-nkBT, 27rfcsr]; temperature T = 4.2 K. 



III. RESULTS AND DISCUSSION 

Figure [2] shows the charge density distribution (a)-(c) 
and the magnctosubband structure (d)-(f) of a GNR for 
several gate voltages, eVg/to = 0, 0.1, and 0.5 respec- 
tively in the regime of the lowest LL, \v\ <\ {v ~ nh/eB 
being the filling factor and n is the average charge con- 
centration). An application of the gate voltage induces 
extra charges in the system, Fig. ^ (b),(c), which leads 
to an asymmetry between the electron and hole branches 
in the energy spectrum, sec Fig. [2](e),(f). For small gate 
voltages, see Fig. [5] (e), the lowest unoccupied spin- up 
branch deforms and gets pinned to the Fermi energy. As 
a result, a region with compressible densities is formed 
in the center of the ribbon. We define the compressible 
strip as a region where the dispersion lies within the en- 
ergy window \E—Ep\ < iksTi^ The compressible strips 
form because in the above energy window the states are 
partially filled, i.e. < fpD < 1 and hence the system 
has a metallic character. The formation of compressible 
strips in graphene has been recently discussed ir>iiiiS. 
Further increase of the gate voltage leads to increase of 
the total electron density, which, in turn, causes even 
stronger pinning of the lowest LL to the Fermi energy 
such that the compressible strip gets extended almost 
over the entire ribbon, see Fig. [2](f). Figures [2] (a)-(c) 
show a density distribution for spin-up and spin-down 



electrons, %(?;), n^{y) and the total density distribution, 
n(i/) = n-|-(y) + n^(y). Due to the strong Zeeman splitting 
the electron branch is predominantly spin-up polarized, 
whereas the hole branch is spin-down polarized. The 
total charge distribution shows a density enhancement 
towards ribbon's edges. This feature is related to the ef- 
fect of the electrostatic Coulomb repulsion in a structure 
with a hard-wall confinement and was discussed (for the 
spinless case) iniiii^. 

Compressible strips in GNRs are characterized by 
the existence of countcrpropagating states at the same 
cdgoi^. This can play an important role in understanding 
of the magnetotransport of interacting electrons in the 
vicinity of the Dirac point. In one-electron picture the 
number of propagating states in the region — 1 < < 1 
always equals to two with electrons of different spins 
propagating near different boundaries, see Fig. [5] (d). 
(Note that a number of propagating states at a given en- 
ergy is given by a number of intersections of the line 
E = Ep with the dispersion bands E = E{k), and 
the direction of the propagation is determined by a sign 
of dE/dkj;.) The backscattering is probable only in 
events accompanied by a spin-flipping since the coun- 
tcrpropagating states with the same spin are located 
at the different edges. However in graphene a spin- 
relaxation length is rather large and can reach several 
micrometers^i, therefore backscattering is expected to be 
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rather unimportant. However, the situation is different 
when electron-electron interaction is taken into account. 
In this case there arc a number of additional transport 
channels at the Fermi energy where electrons of the same 
spin propagate in the opposite directions at the same 
boundary, see Fig. [2](f). In this case no spin-flipping is 
apparently needed for the backscattering and the latter 
can be caused by e.g. impurity scattering of electrons of 
the same spin. To summarize, our self-consistent calcula- 
tions demonstrate that, in comparison to the one-electron 
picture, electron-electron interaction leads to the drastic 
changes in the dispersion relation and structure of prop- 
agating states in the regime of the lowest LL such as a 
formation of the compressible strip and opening of addi- 
tional conductive channels in the middle of the ribbon. 

We proceed to the discussion of the conductance of 
GNRs in the regime of the lowest LL taking into ac- 
count the electron interaction and spin effects. We will 
study the effect of different types of disorder on GNR's 
conductance, focussing on the robustness of respectively 
edge and bulk state transmission. Let us start with the 
case of an ideal GNR (i.e. without any disorder) whose 
conductance is shown in Fig. [3] as a red (solid) curve. In 
the range < < 0.3 there are only two states in the 
energy interval A^nksT and the transport is determined 
by the edge states only, see Fig. [5{e) (note that in this 
range of v the states in the compressible regions practi- 
cally do not contribute to the conductance). As a result, 
the conductance of the ideal ribbon equals to 2Go. Fur- 
ther increase of the filling factor leads to the enhancement 
of the conductance due to the additional states forming 
the compressible strip, see Fig. [2jf). The conductance of 
the ideal ribbon is therefore not quantized and develops a 
bump-like structure with the increase of the filling factor 
above v > 0.3. The suppression of conductance quanti- 
zation in an ideal GNR in a perpendicular magnetic field 
for v > \ was discussed in details ini^. 

Let us now study the effect of impurities on the con- 
ductance, see Fig. [3] (a). Recent two-terminal measure- 
ments show that the quantized conductance plateaux are 
not observed even in high magnetic field^^^— . In particu- 
lar, this behavior is attributed to the effect of impurities 
which can crosslink the chiral currents flowing at opposite 
edgesi^ We utilize the model of short-range impurities 
which are uniformly distributed over the GNR. They are 
modeled by adding to the self-consistent potential a term 
which is randomly chosen in the energy interval [— 15, (5], 
where 5 = 4<o is an impurity strength. The concentration 
of impurities is chosen in the range (0.01 — 0.1)nm~^. As 
expected, the states forming the compressible strip are 
the most sensitive to disorder. They propagate in the 
center of the ribbon and have a significant spatial exten- 
sion in the transverse direction. When the concentration 
of impurities is moderate (green dotted and blue solid 
lines) , the conductance decreases but remains larger than 
2Go since there are at least two edge states which trans- 
mit the current. The edge states are very robust and can 
be suppressed only if the impurity concentration greatly 
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FIG. 3: (Color online). Conductance of the GNR as a func- 
tion of the filling factor v calculated in the presence of differ- 
ent types of disorder: (a) short-range impurities with strength 
5 = 4to, (b) fluctuation of the normal to the surface compo- 
nent of the magnetic field due to the warping; (c) spin-flip 
effect. The red (solid) line on all plots corresponds to an 
ideal ribbon (i.e. without disorder). The arrows (e) and (f) 
indicate the filling factors for which the corresponding band 
structures are plotted on Fig.(l2]| (e) and (f), respectively. The 
inset in (b) illustrates the landscape of a corrugated GNR. A 
typical wavelength of the ripples is about 8 nm, while the 
height fluctuates within 1 nm. T — 4.2 K. 



exceeds typical experimental values (green dashed line). 

Another source of disorder can arise from a surface 
warping which is an inherent property of graphene^^. 
Even though the magnetic field is uniform, a normal to 
the surface component of the magnetic field fluctuates 
due to the corrugated geometry of graphene. It results in 
a spatial-dependent hopping integrals. It has been shown 
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that these spatial correlations can affect the quantum 
Hall (QH) transition at -E = leading to anomalously 
abrupt behaviour— or result in a field-driven topological 
transition from the QH-metal state to the QH-insulator 
state in the vicinity of ~ Oi^ In order to take this effect 
into account we first generate a corrugated surface (see 
inset in Fig. [3](b)) using a method described in Refs^ii^. 
Then graphcnc lattice is projected on the surface and a 
hopping integrals between two adjacent carbon atoms are 
calculated by integration of the vector potential along a 
line linking these atoms according to Eq. ([2|). Wc find 
that the conductance of the warped GNR is practically 
the same as that one of the ideal ribbon, see Fig. [3] (b). 
We therefore conclude that for realistic nanoribbons the 
spatial modulation of the magnetic field due to the warp- 
ing has a negligible effect on the conductance. 

Finally wc investigate the effect of spin-flipping on 
the magnetotransport of GNRs. It is modeled by in- 
troducing in the Hamiltonian an additional term Hsf — 

a J2r ('^lif'^ri + o-^iO-r-t^ , whcrc a describes the local 
spin-flip rate. This term will in generally admix the 
counter-propagating states of different spins leading to 
backscattering. In our calculations the value of a is var- 
ied in the range (0.0— 0.3)tcr^. The conductance of GNRs 
in the presence of spin- flipping is shown in Fig. ^c). 
When a = O.lto (dotted line) the states constituting the 
compressible strips get scattered and the total conduc- 
tance is only due to the edge states. As a increased up to 
0.3to (dashed line) the conductance decreases below the 
value of 2Gq due to stronger scattering but still remains 



mostly dominated by the edge states. Backscattering of 
spin-polarized states due to spin-flipping was suggested 
as the way to manipulate spin currentai^. However our 
results show that even in the case of strong admixture of 
counterpropagating states, the edge state transmission 
remains practically unaffected. 



IV. CONCLUSIONS 

The magnetotransport properties of GNRs in the 
regime of the lowest LL was studied taking into account 
electron interaction, gate electrostatic and spin effects. 
When the gate voltage is applied extra charges are in- 
duced in the ribbon leading to an asymmetry between 
the electron and hole branches in the energy spectrum. 
The lowest LL gets pinned to the Fermi energy and a 
compressible strip is formed in the middle of the GNR. 
There are two types of current carrying states in the low- 
est LL, the conventional edge states and the bulk states 
in the middle of the ribbon constituting the compressible 
strip. The later arc very sensitive to disorder and get 
scattered even if the concentration of disorder is mod- 
crate. In contrast, the edge states are very robust and 
can not be suppressed even in the presence of a strong 
spin-flipping. 
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